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Abstract 



O 

CN \ Maximal signal and peak of high-frequency relic gravitational waves (GW's), recently expected by quintessen- 



tial inflationary models, may be firmly localized in the GHz region, the energy density of the relic gravitons in 



critical units (i.e., 1iq£Igw) is of the order 10 6 , roughly eight orders of magnitude larger than in ordinary infla- 
tionary models. This is just right best frequency band of the electromagnetic (EM) response to the high-frequency 

O ' 

£T) ' GW's in smaller EM detecting systems. We consider the EM response of a Gaussian beam passing through a static 

o" 

magnetic field to a high-frequency relic GW. It is found that under the synchroresonance condition, the first-order 
i . perturbative EM power fluxes will contain "left circular wave" and "right circular wave" around the symmetrical 
?P" axis of the Gaussian beam, but the perturbative effects produced by the states of + polarization and x polarization 

> : 

of the relic GW have different properties, and the perturbations on behavior are obviously different from that 
of the background EM fields in the local regions. For the high-frequency relic GW with the typical parameters 
v g = 10 10 Hz, h = 1CT 30 in the quintessential inflationary models, the corresponding perturbative photon flux 
passing through the region 10" 2 m 2 would be expected to be 10 3 s _1 . This is largest perturbative photon flux we 
recently analyzed and estimated using the typical laboratory parameters. In addition, we also discuss geometrical 
phase shift generated by the high-frequency relic GW in the Gaussian beam and estimate possible physical effects. 
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I. INTRODUCTION 



Relic GW's are very important sources of information on the very early universe, physical behavior 
of the relic GW's expresses the states and evolution of the early universe. Whether direct detection or 
indirect tests to the relic GW's, both of them might provide the new ways to observe our universe. On 
the other hand, the expected properties of the relic GW's, such as amplitudes, polarizations, frequency 
band, energy densities and spectra, et al., are dependent on the concrete universe models. Thus the 
expected features of the relic GW's and the concrete universe models have a closed relation. In recent 
years, the quintessential inflationary models have been much discussed and some astrophysical 

and cosmological observations seem to indicate that AO, El the quintessential inflationary models are 
explicit observationally acceptable. One important expectation III |2J] of the models is that maximal 
signal and peak of the relic GW's are firmly localized in the GHz region, the corresponding energy 
density of the relic gravitons is almost eight orders of magnitude larger than in ordinary inflationary 
models, and the dimensionless amplitude of the relic GW's in the region can roughly reach up to 10~ 30 
[ 1]. This is about five orders of magnitude more than that of the standing GW discussed in Refs.lQ]- 
|2D. Moreover, because the resonant frequencies between the GW's and EM fields in some smaller EM 
detecting systems (e. g., microwave cavities, strong EM wave beams, and so on) [7-12], are just right 
distributed in the GHz region, the results offered new hopes for the EM detection of the GW's (Note 
that the frequency band detected by the VIRGO[13], LIGO (Laser Interferometer Gravitational Wave 
Observatory)[14, 15], and LISA (Laser Interferometer Space Antenna)[16], are often distributed in the 
region of 1CT 4 — 10 4 Hz (this is also most promising detecting frequency band for the usual astronomical 
GW's), thus the EM response to the high-frequency relic GW's might provide a new detecting window 
in the GHz band). 

In this paper, we shall study the EM response to the high-frequency relic GW by a Gaussian beam 
propagating through a static magnetic field. Here we consider it since there are following reasons: 
(1) Unlike the usual EM response to the GW by an ideal plane EM wave[17], the Gaussian beam is a 
realized EM wave beam satisfying physical boundary conditions, and because of the special property of 
the Gaussian function of the beam, the resonant response of the Gaussian beam to the high-frequency 
GW's has better space accumulation effect (see figure 5) than that of the plane EM wave (see fig- 
ure 7 in Ref.[8]). (2) In recent years, strong and ultra-strong lasers and microwave beams are being 
generated[ 18-21] under the laboratory conditions, and many of the beams are often expressed as the 
Gaussian-type or the quasi-Gaussian-type distribution, and they usually have good monochromaticity in 
the GHz region. (3) The EM response in the GHz band means that the dimensions of the EM system 



may reduced to the typical laboratory size (e.g. magnitude of meter), thus the requirements of other 
parameters can be further relaxed. (4) Unlike the cavity electrodynamical response to the GW's (in 
general, the detecting cavities are closed systems for the normal EM modes stored inside the cavities), 
the Gaussian beam propagating through a static EM field is an open system. In this case the EM pertur- 
bations might have a more direct displaying effect, although they have no energy accumulation effect in 
the cavity electrodynamical response. Therefore, the EM response of the Gaussian beams to the GW's 
and the EM detection of the microwave cavities to the GW's have a very strong complementarity each 
other, this is also one of motivations for this investigation. 

The basic plan of this paper is the following. In Sec. II we shall present usual form of the Gaussian 
beam in flat space-time.In Sec.III we shall consider the EM response of the Gaussian beam passing 
through a static magnetic field to the high-frequency relic GW's in the quintessential inflationary models. 
It includes the perturbation solutions of the electrodynamical equations in curved space-time; the first- 
order perturbative EM power fluxes (or in quantum language-perturbative photon fluxes). Moreover, we 
shall give numerical estimations on our results. In Sec.IV we will discuss possible geometrical phase 
shift produced by the high-frequency relic GW's. Our conclusions will be summarized in Sec.V. 

H. THE GAUSSIAN BEAM IN FLAT SPACE-TIME 

It is well known that in flat space-time (i.e., when the GW's are absent) usual form of the fundamental 
Gaussian beam is [22] 

ijjQ 7"^ Z k IT 

^ = + j^jjf exp (~w^ exp ^ keZ ~ ^ ~ tanl j + 2k + (1) 

where r 2 = x 2 + y 2 , k e = 2n/X e , f = 7rW 2 /\ e , W = W [l + (z/f) 2 ]K R = z + f 2 /z. ^ is the 
maximal amplitude of the electric (or magnetic) field of the Gaussian beam, i.e., the amplitude at the 
plane z = 0, W is the minimum spot size, namely, the spot radius at the plane z — 0, and 5 is an 
arbitrary phase factor, ip satisfies the scalar Helmholtz equation 

VV + fceV = 0. (2) 

For the Gaussian beam in the vacuum, we have k 2 = Ld e 2 fi e , where uo e is the angular frequency of the 
Gaussian beam. 



Supposing that the electric field of the Gaussian beam is pointed along the direction of the rr-axis, 
that it is expressed as Eq.(l), and a static magnetic field pointing along the y-axis is localized in the 
region —1/2 < z < 1/2, then we have 
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where the superscript denotes the background EM fields, the notations ~ and A stand for the time- 
dependent and static fields, respectively. Using (we use MKS units) 



B ' = -— VxP, 



(4) 



and Eqs.(l) and (3), we obtain the time-dependent EM field components in the cylindrical polar coordi- 
nates as follows 

Ei 0) = ^cos0, Ef = sin 0, = 0, (5) 
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With the help of Eq.(l) and Eqs.(5)-(8), we can calculate the power flux density of the Gaussian in 
flat space-time. For the high-frequency EM power fluxes, only nonvanishing average values of these 
with respect to time have an observable effect. From Eqs.(l) and (5)-(8), one finds 
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where < S z >, < S r > and < > represent the average values of the axial, radial and tangential EM 
power flux densities, respectively; the angular brackets denote the average values with respect to time. 
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We can see from Eqs.(9)-(11) that < S r > z=0 = < > z=0 = 0, < S z > z=0 = < S z > max and 
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in the region near the minimum spot. Thus, the propagating 
direction of the Gaussian beam is exactly parallel to the z-axis only in the plane z — 0. In the region 

(0) (0) (0) 

of z 7^ 0, because of nonvanishing < S r > and < > and the positive definite property of < S r >, 
Eq.(10), the Gaussian beam will be asymptotically spread as \z\ increases. 

We will show that when the GW is present, the Gaussian beam will be perturbed by the GW. Es- 
pecially, under the synchroresonant condition (i.e., when the frequency uo g /2-n of the GW equals that 



uj e /2n of the Gaussian beam), nonvanishing first-order perturbative EM power fluxes can be produced, 
and they on behavior are obviously different from that of the background EM fields in the local regions. 

III. ELECTROMAGNETIC RESPONSE OF THE GAUSSIAN BEAM TO A HIGH-FREQUENCY 
RELIC GRAVITATIONAL WAVE 

A. The high-frequency relic GW's in the quintessential inflationary models 

For the relic graviton spectrum in the quintessential inflationary models, recent analyses[l, 2] seem to 
indicate that unlike the ordinary inflationary models, the maximum signal and peak are firmly localized 
in the GHz region, the corresponding energy density of the relic gravitons is almost eight orders of 
magnitude larger than that in ordinary inflationary models, and the dimensionless amplitude of the 
relic GW's in the GHz band can roughly reach up to 10 _30 [1]. Thus smaller EM detecting systems 
(not necessarily interferometers) may be suitable for the detection purposes of the high-frequency relic 
GW's. 

In Fig.l one illustrates the relic graviton logarithmic energy spectra expected by the typical ordinary 
inflationary models(curve I) and by the quintessential inflationary models (curve II), respectively. Re- 
gions (1), (2), (3), (4)-l, (4)-2, (5) and (6) represent the detecting frequency bands for the LISA[16], 
LIGO[14, 15], resonant-mass detectors[23-25], the superconducting microwave cavities[10, 11, 38-40], 
the Gaussian beams tuned in on GHz frequency band and Mini-ASTROD (Mini-Astrodynamical Space 
Test of Relativity using Optical Devices) [26], respectively. Present operating-mass detectors include 
ALLEGEO[23], EXPLORER[24], AURIGA and NAUTILUS[25]. These detectors are often operating 
at KHz frequency band. For example, the operating frequency of the cryogenic resonant-mass detector 
EXPLORER is 923 Hz[24]. The sensitivity of them would be expected to be Sh ~ 10 -19 - 1(T 22 
roughly for the GW's of KHz band. 

Giovannini[l] analyzed the relic GW's in the quintessential inflationary models within a framework 
of quantum theory. In this framework the Fourier expansion of the field operators of the relic GW's can 
be written as 

M^V) = J d%[^{k g ,r)) e ik °- s + p,l(k g ,rj) e~^% 

Mw) = j^y^ J d%[jl®(k g ,r)) e ik °- s + jll(k g ,r)) e~ ik °% (12) 

where 
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FIG. 1 : Curves I and II represent the relic graviton logarithmic energy spectra (in critical units) expected by the usual 
ordinary inflationary models and by the quintessential inflationary models, respectively. The curves quote from Ref.[l]. Here 
we illustrate roughly the distribution for some detecting-frequency bands. The region (1) expresses the detecting frequency 
band of LISA (<~ 1CP 4 — 1Hz); the region (2) is one of LIGO (~ 1 — 10 A Hz); the region (3) is one of the resonant- 
mass detectors (~ 10 3 iiz); the region (4)-l is one of the difference-frequency resonant response of the microwave cavities 
(~ 10 3 — 10 4 H z), the region (4)-2 (~ 10 8 — 10 9 Hz) is one of the fundamental resonant response of the microwave cavities; 
the region (5) is one of the Gaussian beams (~ 10 9 - 10 n Hz) and region (6) is one of Mini-ASTROD (~ 1CT 6 - 10~ 3 Hz). 
It can be seen that the region (4)-2 coincides partly with the maximal signal and peak of the relic graviton energy spectra 
expected by the quintessential inflationary models, the detecting frequency band (5) of the Gaussian beam tuned in on GHz 
frequency region can be almost completely localized in this peak value region. For the relevant background in figure 1, we 
presented a brief introduction in Appendixes (Appendix A: The cavity electromagnetic response to the GW's. Appendix B: 
The dimensionless amplitude h and the power spectrum Sy, of the relic GW's. Appendix C: ASTROD. Appendix D: Noise 
problems). 
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© and ® denote + and x polarization of the GW's, respectively, 77 is the conformal time, , and 
a®, a® represent the creation and annihilation operators of the above two physical polarization states. 
From the viewpoint of observation, the classical picture of Eq.(12) corresponds to the amplitudes of the 
relic GW's with the states of the two polarization, namely, 

M^) = j^j^ J d% K{k g ,r,) e ik ~°- s , 

M^) = (2^372 / d3k 9 h ®( k 9>v) e ikg - S , (14) 

where the integration spreads over the whole frequency band of the relic GW's. However, if we hope 
to realize the resonant response of a monochromatic Gaussian beam to the relic GW's, it is necessary 
to let the frequency of the Gaussian beam equals a certain frequency in the peak value region of the 
relic GW's (i.e., uo e = uj g ). Fortunately, most of strong microwave beams (including the Gaussian 
beam) generated by present technology are monochromatic or quasi-monochromatic, thus resonateable 
relic GW to the Gaussian beam will be only one monochromatic component of satisfying condition 
ujg = uo e in the relic GW frequency band. In this case the corresponding treatment can be greatly 
simplified without excluding the essential physical features. Of course, for the resonant response to the 
relic GW's on the earth, we should use interval of laboratory time (i.e., c dt = a(rj) drj) and laboratory 
frequency [27]. Consequently, a monochromatic circular polarized plan relic GW propagating along the 
z-axis can be written as 



h e = h xx = -h yy = A 9 exp(ik a x a ) = A^exp [i(k g z - u g t)], 

h® = h xy = h yx = iA® exp(ik a x a ) = iA® exp [i(k g z - oj g t)\. (15) 

This is just usual form of the GW in the TT gauge. Eq.(15) can be viewed as the classical approximation 
of the Eqs.(12) and (13) under the monochromatic wave condition. 

Since the relic GW's in the peak value region of the quintessential inflationary models have very high 
frequency (the GHz frequency band), future measurements through the microwave cavities or the strong 
Gaussian beams may be useful. Especially, for the high-frequency relic GW's of 10 9 Hz < v g < 10 11 
Hz (these are just best peak frequency band expected by the quintessential inflationary models), using 
the EM response of the Gaussian beam would be more suitable (see Fig.l). 



B. The electromagnetic system in the high-frequency relic gravitational wave field 



From Eq.(15), the nonvanishing components of metric tensor in Cartesian coordinates are given by 

fi'oo — ~ 1? #11 — 9xx = 1 + h xx = 1 + h®, (722 = 9yy = 1 + h yy = 1 — h®, 

9i2 = 921 = 9xy = 9yx = h xy = h®, #33 = g zz = 1, 
g 00 = -1, t? 11 =gT ! = l-h xx = l- K, g 22 = gyy = 1 - h™ = 1 + /i e , 



/ 12 = g 21 = g xy = g yx = -h xy = -h*, g 33 = g zz = 1. 



(16) 



With the help of Eqs.(3)-(8) and considering perturbation produced be the weak GW field expressed as 
Eqs.(15) and (16), the components of the EM field tensor in the Cartesian coordinates can be written as 
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where F$ and represent the background static magnetic field By'' and the background EM wave 
field (the Gaussian beam), respectively, F$ is the first-order perturbation to the background EM fields 
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in the presence of the GW. For the nonvanishing F$ and F$ , we have 



F, 
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The EM response to the GW can be described by Maxwell equations in curved space-time, i.e. 
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where J M indicates the four-dimensional electric current density. For the EM response in the vacuum, 
because it has neither the real four-dimensional electric current nor the equivalent electric current caused 
by the energy dissipation, such as Ohmic losses in the cavity electrodynamical response or the dielectric 
losses[10], so that = in Eq.(18). 

Unlike the interaction of a plane EM wave with a plane GW (according to the Einstein-Maxwell 
equations under the condition of the weak gravitational field, if both the plane GW and plane EM wave 
have the same propagating direction, then the perturbation of GW to the EM wave vanishes[17])., the 
electric and magnetic fields of the Gaussian beam are nonsymmetric (see Eqs.(l) and (5)-(8)). In this 
case, using Eqs.(5)-(8) and (18)-(19), it can be shown that even for the plane GW propagating along 
the positive ^-direction, it can produce nonvanishing perturbative effect to the Gaussian beam. In order 
to find the concrete form of the perturbation produced by the direct interaction of the GW with the 
Gaussian beam, it is necessary to solve Eqs.(18) and (19) by substituting Eqs.(l), (5)-(8), (15), (16) and 
(17) into them, which is often quite difficult. However, as shown in Refs.[8, 17, 28, 29], the orders of 
amplitudes of the first-order perturbative electric and magnetic fields produced by the direct interaction 
of the GW with the EM wave (e. g., plane wave or the Gaussian beam) are approximately hB^c and 
hB(°\ respectively, while that generated by the direct interaction of the GW with the static magnetic 
field are approximately hB^c and hB^ respectively. Thus corresponding amplitude ratio is about 
hB^/hB^. In our case, we have chosen B^ ~ 1(T 3 T, ~ 10T, i.e., their ratio is only 1(T 4 
roughly. Therefore, the former can be neglected. In other words, the contribution of the Gaussian 
beam is mainly expressed as the coherent synchroresonance (i.e., uo e = oo g ) of it with the first-order 
perturbation F$ generated by the direct interaction of the GW with the static field By°\ In this case, 
solving the process of Eqs.(18) and (19) can be greatly simplified, i.e., the static magnetic field B^ = 
By ^ can be seen as the unique background EM field in Eqs.(18) and (19). Under these circumstances, in 
the first we can solve Eqs.(18) and (19) in the region II (-1/2 <z<l/2, B<® = Bf ] ) to find the first- 
order perturbation solutions; and second, using the boundary conditions one can obtain the first-order 
perturbation solutions in the region I (z < -1/2, B^ = 0) and the region III (z > 1/2, B^ = 0). 

Introducing Eqs.(15)-(17) into Eqs.(18)-(19), neglecting high-order infinite small quantities and 
perturbative effect produced by the direct interaction of the GW with Gaussian beam, Eqs.(18) and (19) 
are reduced to 
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and 

E% = E$ = 0, B&> = BiJ = 0, (22) 

where the commas denote partial derivatives. Using Eq.(15), Eqs.(20)-(22) can also be expressed as 
following inhomogeneous hyperbolic-type equations, respectively: 

= Ei% - ±E% = -A e B^k 2 g c exp - u g t)}, 
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= B^ z - Is« = i^B^ exp [i( V - (24) 

= = 0, (25) 

where □ indicates the d'Alembertian. Obviously, every solution of Eq.(20) must satisfy Eq.(23), and 
every solution of Eq.(21) must satisfy Eq.(24), and it is easily seen from Eqs.(22) and (25), that a 
physically reasonable solution of them would be only 

m = BP = 0. (26) 



The general solutions of Eqs.(20)-(24) in the region II (—1/2 < z < 1/2) are given by 

E^ = l-A^B^kg c(z + 1/2) exp [i(k g z - u g t)] + bi exp [i(k g z - u g t)] + c x exp [i(k g z + w g t)], 

B<f> = l -A e B^k g (z + 1/2) exp [i(k g z - u g t)\ + b 2 exp [i(k g z - u g t)\ + c 2 exp [i(k g z + u g t)\, (27) 
E^ = —AtgBjjpkg c(z + 1/2) exp [i(k g z - uj g t)\ + ib 3 exp [i(k g z - uj g t)\ + ic 3 exp [i(k g z + w s t)], 

S^. 1 ^ = ^A®B^k g (z + //2) exp [i(fc g ^ - u g t)\ + ib 4 exp [i(k g z - u g t)\ + ic 4 exp [i(k g z + u g t)]. (28) 

The solutions (26)-(28) show that perturbative EM waves produced by the plane GW satisfing TT 
gauge, must be the transverse waves, and the background static magnetic field is perpendicular to the 
propagating direction of the GW. The constants in Eqs.(27) and (28) satisfy 

6i - cb 2 = --AqBWc, ci + cc 2 = 0, 

b 3 + cb 4 = - l -A^c, c 3 - cc 4 = 0, (29) 

their concrete forms will be defined by the physical requirements and boundary conditions. 

The solutions (26)-(28) have the similar features as that found by a number of authors [28, 29] 
previously, but as we shall show that unlike the previous works, our EM system will be a possible 
scheme to display the first-order EM perturbations produced by the high-frequency relic GW's (GHz 
region) inside of the typical laboratory size (not necessarily interferometers or EM cavities with giant 
dimensions), a particularly interesting feature of the first-order perturbation will be the perturbative 
effect in the some special directions and in the some special regions. 

It should be pointed out that in curved space-time, only local measurements made by an observer 
travelling in his world-line have definite observable meaning. These observable quantities are just the 
projections of the physical quantities as the tensor on tetrads of the observer's world-line. The tetrads 

1 o 



consist of three space-like mutually orthogonal vectors and a time-like vector directed along the four- 
velocity of the observer, the latter is perpendicular to the preceding ones. We indicate these with r^, 
where the index in brackets numbers the vectors and the other refers to the components of the tetrads 
in the chosen coordinates. Consequently, the quantities F^ a ^ measured by the observer are the tetrad 
components of the EM field tensor, that is 

= FprTfaTfa. (30) 

Obviously, for our EM system, the observer should be laid at rest state to the static magnetic field, i.e., 
only the zeroth-component of his four-velocity nonvanish. Thus, the tetrad has the form 

7f 0) = (r ( ° 0) , 0, 0, 0). (31) 

Using Eq.(16) and the orthonormality of the tetrads d^uT^T^ = rf a p, neglecting high-order infi- 
nite small quantities, it is always possible to get 

rg } = (1, 0, 0, 0), 

= (0, 1 - l -h m 0, 0), 

T (2) = (°> ~ h ®i 1 + \ h ®, 0), 

rg } = (0, 0, 0, 1). (32) 

Eq.(32) indicates that the zeroth and third components of the tetrads coincide completely with the 
time- and z-axes in the chosen coordinates, respectively. Furthermore, has only the projection on 
the x-axis, thus points actually at the rr-axis. This means that the azimuth <fi in the tetrads and that in 
the chosen coordinates are the same. While the deviation of from the |/-axis is only an order of h®. 

With the help of Eqs.(3),(4),(27),(28) and (32), neglecting high-order infinite small quantities, we 
obtain 
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As we have pointed out above, here \hF$\ ~ AB^ , \F$\ ~ AB y 0) , and in our case, 

B (o) 

10T, ~ 1(T 3 T, thus we have (/iF^Vl^l ~ B^/B y 0) « 1(T 4 , so that /iF^? terms in 
Eq.(33) can be neglected again. In this case Eq.(33) can be further reduced to (in the region II: — Z/2 < 

z < 1/2) 
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+ c 2 exp [i(k g z + uj g t)], (34) 



= ~\A®BfX c (z + Z/2) exp [<( ^ - Ug t)] + ife exp [i( ^ - Ml + <P3 exp [i(k g z + 
S (x) = + ^ 0J B(°) 

= \A®Bfkg{z + Z/2) exp [i( - + exp [i(k a z - ^)] + ^ exp [i(k g z + Wfl t)] 

exp - w fl f)], (35) 

and 

E (z) = 0, 

^) = ^ 0) = -^ • (36) 
C. The particular solutions satisfying boundary conditions 

In fact, the perturbative parts in Eqs.(34)-(36) are the general solutions of Eqs.(20)-(25) in the region 
II (—Z/2 < z < 1/2). We shall define the constants in Eqs.(34) and (35) to give the corresponding 
particular solutions satisfying the boundary conditions. 



Clearly, the perturbative EM fields in the regions I, II and III must satisfy the boundary conditions 
(the continuity conditions): 

(■^(£oj)*=-i/2 = (F$)ii)*=-i/2 ' (-^S)j/)«=»/2 = (^>om)*=i/2 • ( 37 ) 
If one chooses the real part of the pure perturbative fields in Eqs.(34) and (35), then we have 

MH = —A^B^kg c(z + 1/2) sin (k g z - u g t) + h cos (k g z - u g t) + c x cos (k g z + u g t), 
B (l) = ~7>A(BBf ) k g (z + 1/2) sin (k g z - u g t) + b 2 cos (k g z - u g t) + c 2 cos (k g z + u; fl £), (38) 

F tv) = ~\A®B i y ) k g c(z + 1/2) cos - uj g t) - b 3 sin (& 9 2: - uj g t) - c 3 sin (k g z + w g t), 

—A®B^ sin (k g z — uj g t). (39) 

A physically reasonable requirement is that there is no the perturbative EM wave propagating along 
the negative ^-direction in the region 111(2 > 1/2). Here we shall consider one more simple case, i.e., 
the perturbative EM wave in the negative ^-direction is also absent in the region l(z < —1/2). In order 
to satisfy the boundary conditions, Eq.(37), and the above requirement, from Eqs.(29) and (37)-(39), 
one finds 

b^—A^c, b 2 = \A^Bf\ 6 3 = ^fc, h = - 3 -A^\ 



(40) 

and 

I = n\ g (n is integer). (41) 

In this case we have the perturbative EM fields in the above three regions as follows: 
(a) Region I (z < -1/2, B<® = 0) 

_ c pW - e {1) - cF {1) - 

^(x) ~~ ^(01) ~~ ^(y) ~ Cr {02) ~ U > 

5(1) _ 5(1) _ 5(1) _ 5(1) _ n , AO s 
B ( X ) ~ F (32) ~ B (y) ~ F (13) ~ - ( 42 ) 



(b) Region II (-1/2 <z< 1/2, B^ = B y 0) ) 



1 1 

E^ = --A^Bfkg c (z + 1/2) sin (k g z - uo g t) - -A e B y 0) c sin (k g z) sin (u g t), 

--A^B^k g (z + l/2) sin (k t) + -A Q B<® sin (k g z) sin (UJgt), 



(1) 

(y) 



(43) 



E$ = --A^B^kg c(z + 1/2) cos (k g z - u g t) - ~A®Bf ] c sin (k g z) cos (u g t), 

B{1\ = l -A®B^k g {z + 1/2) cos (k g z - u g t) - \a®B^ sin (k g z) cos (u g t), (44) 

(c) Region III (1/2 < z < l , B^ = 0) 

E$ = ~A 9 B^k g clan(k g z-u g t), 

Sgj = -^Bffc s Z sin (fc fl * - (45) 

£7$ = --A^B^kgdcoaikgZ-Ugt), 

= Z cos - Wfl *)> (46) 

where / is the size of the effective region in which the second-order perturbative EM power fluxes, such 
as ^(E^B^), ^(E^B^), keep a plane wave form. Notice that the power fluxes in the region II 
contain the parts with space accumulation effect, i.e., they depend upon the square of the interaction 
dimension. This is because the GW's and EM waves have same velocity, so that the two waves can 
generate an optimum coherent effect in the propagating direction. It is easily to show that if we choose 
the imaginary part of the pure perturbative fields in Eqs.(34) and (35), we can obtain the similar results, 
but Eq.(41) will be replaced by 

Z = (2n+1) — (n is integer). (47) 
2 



Logi and Mickelson [29] used Feynmann perturbation techniques to analyze the perturbative EM waves 
(photon fluxes) produced by a weak GW (gravitons) passing through a static magnetic (or an electro- 
static) field, and found that perturbative EM waves (photon fluxes) propagate only in the same and in 
the opposite propagating directions of the GW(gravitons), the latter is weaker than the former, or the 
latter is absent. Obviously, our results and calculation by Logi et.al are self-consistent. However, due 

1 c 



to the weakness of the interaction of the GW's(gravitons) with the EM fields(photons), we shall focus 
our attention on the first-order perturbative power fluxes produced by the coherent synchroresonance 
of the above perturbative EM fields with the background Gaussian beam, not only the second-order 
perturbative EM power fluxes themselves, such as ^(E^B^]) and ^(E^B^). 



D. The first-order perturbative electromagnetic power fluxes 



The generic expression of the energy-momentum tensor of the EM fields in the GW fields is given 



by 



Ho V 4 / 



(48) 



Due to F^ = F$ + F$ and 
disintegrated into 



=.(i) 



« 



7 (0) 



for the nonvanishing F$ and F$ , T^ v can be 



(0) (1) (2) 

pflU pfMU _|_ pflV _|_ rpfiV 



(49) 



(0) (1) (2) 

where T^ v is the energy-momentum tensor of the background EM fields, and T^ v and T^ v represent 

(0) 

first- and second-order perturbations to T^ v in the presence of the GW, respectively. 
Using Eq.(16), T M ", T M " and can be written as 



(°) i / 1 \ 

no v a ^ a " )' 



(50) 



(i) 
pw 



1 

Ho 



- (p^F ua{1) + FWF" a (°)j + It/ 4 " ^F^F q/3(0) + F^F C 



■m\ _ -h^F {0] F am 
) 4 a/3 



(51) 



(2) i 



(52) 



(0) (1) 



(2) 



For the nonvanishing T , T and T , we have 



(0) (1) (2) 

pfll/ I I p/IU I I pfll/ I 



(53) 



I "7 



Therefore, for the effect of the GW, we are interested in but not in T 1 *" and T M \ Nevertheless, 

(2) o 

it can be shown from Eqs.(l), (3), (5)-(8),(15), (42)-(46) and (52), that the average value of T°° with 
respect to time is always positive. Thus it expresses essentially the net increasing quantity of the energy 
density of the EM fields. Especially, under the resonant conditions, it should correspond to the resonant 

(2) (1) 

graviton-photon conversion at a quantum level [29, 30, 31]. But because of | T^ u \ « \ T^ u | for 

(1) (2) 

the nonvanishing T^ v and , the second-order perturbations are often far below the requirements 
of the observable effect. In this case it has only theoretical interest. However, for some astrophysical 
situations, it is possible to cause observable effects, because where very large EM fields and very strong 
GW's often occur simultaneously and these fields extend over a very large area [32, 33]. 
By using Eqs.(l), (3), (5)-(8), (15), (42)-(46) and (51). we obtain 



S> = = i (mm) = -i (mm i («) «•*, m 



£ = c# = - J. ( £ <« B <°>) = _ JL (^» B <».) cos ^ _ J. («) sm * (55) 



Si=c^ = l(Eg»as») + l(E£)agj). (56, 



(1) (1) (1) 

where ,5"^ and S* 2 represent the first-order radial, tangential and axial perturbative power flux densi- 
ties, respectively. As we have shown above, that for the high-frequency perturbative power fluxes, only 
nonvanishing average values of them with respect to time have observable effect. It is easily seen from 
Eqs.(l), (5)-(8) and (42)-(46), that average values of the perturbative power flux densities. Eqs.(54)- 

(56),vanish in whole frequency rang of uo e ^ uo g . In other words, only under the condition of u e = uo g 

(i) (i) (i) 

(synchroresonance), S r , S 4 ' and S z have nonvanishing average values with respect to time. 

(i) 

In the following we study only the tangential average power flux density < > £i)e=£i , . Introducing 
Eqs.(l), (8), (42)-(46) into (55), and setting 5 = n/2 in Eq.(l) (it is always possible), we have 

(i) (i) (i) 

< S 4 > We =w 9 =< Sq >uj e=U j g + < >ui fi =LL)gi (57) 

where 

(i) i 
<^>^=--<^ ) ) 4°) ) >cos0, 



v ' 1 

< S* >^ =u) = < mm > sin 0. (58) 

(1) (1) 

< S® >ui e =u g and < > u , e .=w g represent the average values of the first-order tangential perturbative 
power flux densities generated by the states of + polarization and x polarization of the GW, Eq.(15), 
respectively. Using Eqs.(l), (8), (42)-(46), (57), (58) and the boundary conditions [see, Eqs. (37) and 
(41)], one finds 

(a) Region I (z < -1/2), 

(i) (i) 

< >U) e =Ul g = < S® >U) e =U)g= 0. (59) 



(b) Region II (-1/2 <z< 1/2), 



(i) 

<st> 



^VggVoV (z + 1/2) 
Vo [l + (z/f)T (z + P/z) 
AgggV (g + Z/2) 

n2i3/2 



,2 Aj f z 
cos I tan -1 — — 



/ 2/2 



^ k ir^ 
sin I tan -1 — — 



4^ [l + (z//) 2 ]*'' V / 2i? 



8/io [l + (^//) 2 ] 1/2 {z + PI 



sin (A; 9 ;z) cos ^tan 1 — 



27? *** 



4^k g W* [I + (z/f) 2 ] 



21 3/2 



tan" 1 - - -|^- - k g z\ j exp(-— ) sin (20) , (60) 



(i) 

< St, > L 



A®Bj 0) jj k g r (z + 1/2) 
4/io [l + (^//) 2 ] 1/2 (^ + / 2 A) 



sin 



A; 3 r 2 
~2~R 



A fj; u V o r (z + 1/2) (k g r 2 



2^W 2 [1 + (z/f) 2 ] 



tan^ 



tan — 



./ 



■ sin (k g z) sin ^A: 9 z — tan 1 j + 3 



4^ [l + (^//) 2 ] 1/2 (^ + / 2 A) 



2i? 



2pt k g W 2 [I + (z//) 2 ] 



21 3/2 



£; 9 z - tan" 1 y + J j exp(-— ) sin 2 0. (61) 



(c) Region III (7/2 < z < l Q ), 



(i) 



A @ B^ Q k g lr 



8fM> [l + (z/f) 2 ] 



21 1/2 



A 9 Bl 0) il, lr 



4^W 2 [1 + (z/f) 2 ] 



21 3/2 



sm 



tan 1 — 



(z + P/z) 
z 

1 



cos tan — 



z 



k„r" 



k g r 2 
~2R 



exp( 



W 2 



f 2R 



)sin (20), 



(62) 



(i) 



4/iQ [i + (*//r] 



21 1/2 



sm 



2^ W 2 [l + (z/f) 2 f 2 



cos 



2# 



(z + f 2 /z) 

2 



-1 -2 

tan — 



./ 



k g r 2 



exp( 



tan^ 



•^ )8in20 - 



(63) 



It is easily shown that the nonvanishing first-order perturbative power flux densities are much greater 
than corresponding second-order perturbative power flux densities. The quantum picture of this process 
can be described as the interaction of the photons with the gravitons in a background of virtual photons 
(or virtual gravitons) as a catalyst [29, 34], which can greatly increase the interacting cross section 
between the photons and gravitons. In other words the interaction may effectively change the physical 
behavior (e.g., propagating direction, distribution, polarization and phase) of the photons in the local 
regions, even if the net increasing quantities of the photon number (the EM energy) of the entire EM 

system approach zero, such properties may be very useful to display very weak signals of the GW's. 

(i) 

Eqs. (60)-(63) show that because there are nonvanishing (which depends on the + polarization 

(i) 

state of the GW) and S® (which depends on the x polarization state of the GW), the first-order tangential 

perturbative power fluxes are expressed as the "left circular wave" and "right circular wave" in the cylin- 

(i) 

drical polar coordinates, which around the symmetrical axis of the Gaussian beam, but < >ui e =ui g 
(i) 

and < >ui e =u g have a different physical behavior. By comparing Eqs.(60)-(63) with Eqs.(9)-(1 1), 

(!) (0) 

we can see that: (a)< > We=Wg and < > have the same angular distribution factor sin(24>) , 

(i) ' (0) 

thus < S$ >uj e =w g will be swamped by the background power flux < >. Namely, in this case 

(i) ' " (i) 

< S® > LUe =ujg has no observable effect, (b) The angular distribution factor of < > We=Wg is a sin 2 <f) 

(0) (i) 

, it is different from that of < >. Therefore,< > We=Wg ,in principle, has an observable effect. In 

(o) (i) (i) 

particular, at the surfaces <fi = 7r/2,37r/2, < S* >= 0, while |< > Ue=Wg \ = |< <S® > We =a; g |max , this 



(0) « 

is satisfactory (although < S r > and < S® >uj e =Lu a have the same angular distribution factor sm 2 0, the 

(0) C 1 ) (0) 

propagating direction of < S r > is perpendicular to that of < 5| > We = Wfl , thus < 5 r > (including 

(0) 

< S z > ) has no any essentially contribution in the pure tangential direction). 

(i) 

Figure 2 gives the distribution of < S§ > We=ulcj at the plane z = 1/2 = |A S (n is a integer) in the 

(i) 

cylindrical polar coordinates, while figure 3 gives the distribution of < > u)(1=UJg at the plane z = in 
the cylindrical polar coordinates. 



" Right circularly " 
propagating region 

0± <S|> 



" Left circularly " 
propagating region 

of <S|> 




Left circularly " 
propagating region 

of <S|> 



" Right circularly " 
propagating region 

of <Sl> 



(i) 

FIG. 2: Distribution of < S@ >u> e =u g at the plane z = 1/2 = ^X g (n is integer) in the cylindrical polar coordinates. It has 
maximum at <f> = 7r/4, 37r/4, 57r/4, and 77r/4, while it vanishes at <f> = 0, 7r/2, 7r and 37r/2, here Z = 0.1m,A 9 = 0.01m and 
r = 0.05to, and the GW propagates along the z-axis. 



(0) 



(0) 



(1) 

<t> 



From Eqs.(10),(ll),(60) and (61), we can also see that < S r >=< S 4, >=< S* > 
plane z — 0, while where 



at the 



(i) 

<st> 



U> e =UJg 



4/ioW^ 2 



exp | ) sin 2 



(64) 



(o) (i) (i) 

In Table I we list the distribution of < > , < > We =nj g and < >^ e =^ s in some typical 

regions. 

Table I, Eqs.(60)-(63) and Figures 2, 3 show that the plane z = and the planes <p = ir/2, 3n/2 

(0) (0) C 1 ) 

are the three most interesting regions. For the former, < >=< S r >=< > Wc =w g = 0, but 




(1) 

FIG. 3: Distribution of < St > u<1 =u>g at me plane z — in the cylindrical polar coordinates. It has maximum at </> = tt/2 

(i) 

and 37r/2, while it vanishes at <f> = 0, tt, Unlike figure 2, < St > We=Wg at the z = is completely "left-hand circular", here 
r = 0.05m, \ g = 0.01m, and the GW propagates along the z-axis. 



(0) (1) (1) 

TABLE I: The distribution of < >,< 5* > We=a , o and < >u e =u B in some typical regions 



Angular distribution 
factor 



The plane 

z = 



The planes 

(f> = 71-/4,5^/4 



The planes 

= 37I-/4, 7tt/4 



The planes 

4> = 7r/2,37r/2 



(0) 

< > 



sin(2(j>) 



(0) (0) 

max ■ 

"Left circularly" 
propagation. 



(0) 

< S"^ >| 



(0) 

< S"> >| max ) 



"Right circularly" 
propagation. 



(i) 



sin(2<p) 



(i) 
(l) 

"Left circularly" 
propagation. 



(i) 

< 5 ffi > We=Wa | : 
(1) 

"Right circularly" 
propagation. 



(i) (i) 
< 8* > We=ulg sin^ |< > Ue=Ug | = 

(l) 

1^ ~^>uj e —ujg | max 

at0 = tt/2, 3tt/2 



There are 

nonvanishing values, 
"Left circularly" 
propagation. 



There are 

nonvanishing values, 
"Left circularly" 
propagation. 



(i) 

(1) 

l< si > 



u) e =u)n I max 



"Left circularly" 
propagation. 



(1) (0) (0) 

nonvanishing < > We =w 9 exists; for the latter, < >=< S@ > u)e=0Jg = 0, but there is non-zero 

(i) 

< > Ue =u g - This means that any nonvanishing tangential EM power flux in such regions will express 
the pure electromagnetic-gravitational perturbation. 



E. Numerical estimations 



If we describe the perturbation in the quantum language (photon flux), the corresponding perturbative 

(i) ~ " (i) 

photon flux caused by < > We =w 9 at the plane <fi = ix/2 {we note that < > We=Wg is the unique 

nonvanishing power flux density passing through the plane) can be given by 

(i) 

n,= % LT J i/2 < S t>^,^/2dzdr, (65) 

(i) { (i) 

where < u% > We=UJg , ^ =7r /2= 

Jo ' J-i/2 < ^® >Lu e =uj g , (j>=ir/2dzdr is the total perturbative power flux 
passing through the plane <fi = n/2, h is the Planck constant. 

In order to give reasonable estimation,we choose the achievable values of the EM parameters in 
the present experiments: (l)V'o = 3 x lO^V.m^ 1 (i. e., ipo/ c — 10~ 3 T), the amplitude of the Gaus- 
sian beam. If the spot radius W of the Gaussian beam is limited in 0.05m, the corresponding power 

w (0) 

can be estimated as P = J Q < S z > z= o 2nrdr 10 5 W (see Eq. (9)), this power is well within the 
current technology condition [21, 35]. For the Gaussian beam with v e = 3 x 10 10 Hz, this is equiv- 
alent to a photon flux n (0) of 3 x 10 28 s _1 roughly. {2)B y 0) = 30T, the strength of the background 
static magnetic field, this is achievable strength of a stationary magnetic field under the present exper- 
imental condition [36]. (3) A® = 10~ 30 , cu g /2n = v g = 3 x 10 10 Hz, these are the typical orders 
expected by the quintessential inflationary models [1]. Substituting Eqs.(61),(63) and the above pa- 
rameters into Eq.(65), and setting W = 0.05m,/ = 0.1m, l = 0.3m,we obtain ps 1.57 x 10 3 s~ 1 
(see Table II). Up to now, this is largest perturbative photon flux in a series of results [7, 8, 9]. Re- 
cently, we analyzed and estimated them under the typical laboratory parameter conditions. If the inte- 
gration region of the radial coordinate r in Eq.(65) is moved to W < r < r (here W = 0.05m, 

r = 0.1m), in the same way, the corresponding perturbative photon flux n'^ can be estimated as 

(i) 

n '<f> = ltJZJ-i/2 <S t>^ a ,^/2dzdr = 0.9Qxl0h- 1 » lO 3 ^ 1 . Although then n' < n , 
it retains basically order of 10 3 s~\ and because the "receiving" plane of the tangential perturbative 
photon flux is already moved to the region outside the spot radius W of the Gaussian beam, it has more 

realistic meaning to distinguish and display the perturbative photon flux. 

(i) (0) 

Figure 4 gives the rating curve between < > CUe =cu g and < S z > at the plane z — 0, here r is the 
radial coordinate. 

Figure 5 gives the rating curve between and the axis coordinate z, and relative parameters are 
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FIG. 4: Rating curve between < > We=Ws and < S z > at the plane z = 0, hers r is the radial coordinate, and <j> = n/2 

(i) (i) 

or 3ir/2 in Eq.(64). It shows that | < > C j c=0 j g \ = \ < > 0Jc=UJg | max at r = Wq/V2. The background axial power 

(0) (0) 

flux density < S z > is typical Gaussian distribution, and | < S z > \ » \ < S@ > We=Wg | ma x- However, because the 

(0) (!) (0) 

propagating direction of < S z > is perpendicular to that of < > We=a)£J , < S z > has no any contribution in the pure 

(i) 

tangential direction. The shaded part expresses < S 1 ^ >„ t=1Js . Here v e = v g = oj g /2ir = 3 x 10 10 Hz, A® = 10~ 30 , 
= 30T, Vo = 3 x ltfV.m- 1 , 1 = 0.1m, W = 0.05m. 



chosen as v e = v g = u 9 /2ti = 3 x 10 l0 Hz, A® = 1(T 30 , = 30T, ip = 3 x 10 5 V.m _1 , / = 0.1m, 
/ = 0.3m and Wo = 0.05m. Figure 5 shows that has a good space accumulation effect as z increases. 



(i) 

In fact, in addition to the third- and fourth-terms in Eq.(61), the rest in the expression of < > 



U) e =U)g 1 



Eq.(61) and Eq.(63), is all slow enough variational function at the z direction. This means that the value 

(i) 

of < > We=Wg is slow variational and keeps its sign invariant in the whole region of the coherent 

(i) 

resonance (here it is about the region of 80A S , namely 0.8m), that they make < > We =uj 3 being all 



"left circularly" propagated from —1/2 to l 

(s _1 ) 
1500 

1250 




z (xlO~Si) 



FIG. 5: Rating curve between and the axial coordinate z, here v e = v g = u) g /2-K = 3 x 10 10 Hz, A® = 10~ 30 , 
By ' 1 = 30T, Vo = 3 x 10 5 V.m -1 , and Wq = 0.05m. It shows that the pure tangential perturbative photon flux passing 
through the plane (f> = tt/2 with 10~ 2 m 2 , would be expected to be 1.57 x 10 3 s~ 1 . 



Table II gives the tangential perturbative photon fluxes and corresponding relevant parameters in 
the three cases. In the first case W 01 = 0.05m, in the second case W 02 = 0.02m, in the third case 
W03 = 0.1m, but the background Gaussian beam keep the same power in such three cases, i. e., 
P pa 10 5 W. Due to W 03 > W i > W 02 , so that ip 03 < ip 01 < ^ 2- Table II shows that the tangential 
perturbation in the Gaussian beam with larger W (i. e., smaller 9) has a better physical effect than 
that in the Gaussian beam with smaller W (i. e., larger 9), here 9 is the spreading angle of the beam, 
e = tarr 1 (-fer) 



nWo 



We emphasize that here oc ip oc v P (see Eqs.(61), (63), and (65), P is the power of the Gaussian 
beam), and at the same time, oc By°\ Therefore, if P is reduced to lO 3 ^, then pa 1.57 x 10 2 s _1 ; 
and even if P is reduced to 10W (this is already a very relaxed requirement), we have still ~ 
1.57 x 10s _1 . Thus, if possible, increasing (in this way the number of the background real photons 
does not change) has better physical effect than increasing P. According to the above discussion, we 
give some values for the power P of the background Gaussian beam and corresponding parameters 



TABLE II: The tangential perturbative photon fluxes and corresponding relevant parameters 



(i) 

A v g (Hz) W {m) 9 <u*>„ e=u)g (W) n^ 1 ) 

(1) 1(T 30 3 x 10 10 0.05 6.36 x 10~ 2 3.11 x 10~ 20 1.57 x 10 3 

(2) 10- 30 3 x 10 10 0.02 15.78 x 10~ 2 1.31 x lO" 20 6.61 x 10 2 

(3) 10~ 30 3 x 10 10 0.10 3.18 x 10~ 2 3.14 x 10~ 20 1.58 x 10 3 



and 71^ (see Table III). 

TABLE III: The power of the background Gaussian beam and corresponding and n' 4 

P(W) n^s- 1 ) n'^s- 1 ) 

10 5 ~ 1.57 x 10 3 ~ 10 3 

10 3 ~ 1.57 x 10 2 ~ 10 2 

10 ~ 1.57 x 10 ~ 10 



In particular, since n'^ indicates the tangential perturbative photon flux passing through the "receiv- 
ing" plane <fi = n/2 (~ 10 _2 m 2 ) outside the spot radius W of the Gaussian beam, the results provided 
a more realistic test scheme. 

IV. GEOMETRICAL PHASE SHIFT PRODUCED BY THE HIGH-FREQUENCY RELIC GRAVITA- 
TIONAL WAVES 

Mitskievich et al. [37] investigated Berry's phase shift of a monochromatic EM wave beam in a 
plane monochromatic GW field, and it is shown that (a) For parallel the propagating GW and EM wave, 



the phase shift is absent, (b) When the both waves are mutually orthogonal, the nonvanishing phase shift 
will be produced, and the phase shift is proportional to the distance propagated by the EM wave (see, 
Eq.(6) in Ref. [37]). i. e., 

Aa = A (A - sin A) sin (A + 5') , (66) 



where 

A = 27CL/\g = ^Ug, (67) 

A is the amplitude of the GW (e. g., A e or L is the distance between the observer and reflecting 
system of the EM wave, and A e = X g , oo e = u g , L is also the interacting dimension of the GW with the 
EM wave. If A >> 1, from Eq.(66), we have 

Aa « AA sin (A + 5') , (68) 



namely, one obtained a linear increase of the phase shift. Eqs.(67) and (68) show that large distance L 
and high-frequency cu g will produce better physical effect than that of the small L and low frequency uj g . 

As we have pointed out that unlike an ideal plane monochromatic EM wave, the Gaussian beam 
is a realized EM wave beam satisfying physical boundary conditions, and as we shall show that for a 
Gaussian beam with a small the spreading angle, the wave beam in the region near the symmetrical 
axis can be approximately seen as a quasi-plane wave. In this case we are possible to estimate the 
geometrical phase shift in the Gaussian beam. 

In order to simplify our analysis, we consider only the real part of Eq.(l). From Eq.(l) we have 

Re{jp) = fiij'i z)cos(k e z — u e t) + f 2 (r, z)sin{k e z — u e t), (69) 



where 
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M'>*) = [l + {z/f fY /2 eM ^ ] f~~2R~ 5 )' (71) 

Eqs.(69)-(71) show that for the Gaussian beam with the small spreading angle, deviation of the propa- 
gation direction from the z-axis in the region near axis would be very small, and for the high-frequency 
band, the functions f\ and / 2 , Eqs.(70) and (71), will be slow variational functions at the z direction. 
In this case the change of the Gaussian beam as space-time mainly depends on the propagation factors 
cos(k e z — uj e t) and sin(k e z — u e t). In this sense it is just characteristic of the plane wave. Therefor, for 
the GW propagating along the y-axis, because its propagating direction is perpendicular to the z-axis, 
it would be to generate the phase shift satisfying approximately Eq.(66) or Eq.(68). Notice that since 
in this case the propagating direction of the GW is parallel with the static magnetic field , whether 
from the classical or the quantum theories of the weak fields, the GW (gravitons) does not produce per- 
turbative EM fields (photon fluxes ) to the static magnetic fields [28, 29]. Thus then the all first-order 
perturbations expressed as Eqs.(54)-(56) vanish. 

It is interesting to compare the geometrical phase shift produced by the high-frequency relic GW of 
v g = 3 x 10 10 Hz and A = 10~ 30 in the Gaussian beam with the phase shift generates by the expected 
astronomical GW of v g = 3 x 10 3 Hz and A = 10~ 22 in the plane monochromatic EM wave (see Table 
IV). With the help of Eqs.(66)-(68), we list the some typical parameters in Table IV. 

We can see from Table IV that scheme (a) has typical parameters of the expected astronomical GW's 
[see analysis in Ref.[37] and Eqs.(60)-(68)], and A = 2.5 x 10 4 , Act = 2.4 x 1(T 18 . This means that in 
this case a GW of A = 10~ 22 and v g = 10 3 Hz can be treated as an effective magnitude of some 1CT 18 , 
but it needs the interacting dimension of cislunar distance (i. e., the reflecting system of the EM wave 
beam is placed on the surface of moon). For the scheme (d), the phase shift may achieve the same order 
of magnitude with the scheme (a), but it needs a interacting dimension of the LISA (~ 10 9 m), and it 
is necessary to constructing a very strong Gaussian beam with the small spreading angle in microwave 
frequency band for the LISA. This seems to be beyond the ability of presently conceived technology. 
Nevertheless, if the amplitude of the high-frequency relic GW ( v g ~ 10 10 Hz) and the amplitude of 
the expected astronomical GW (u g ~ 10 3 Hz) have the same order of magnitude, then the geometrical 
phase shift Act produced by the former will be seven orders larger than that generated by the latter. 



TABLE IV: The geometrical phase shift produced by the GW's in the EM wave beams with the same frequencies 



A u g = u e {Hz) The wave form The interacting A The geometrical 

dimensions L (m) phase shift Aa 

(a) 1CT 22 3 x 10 3 Monochromatic 3.8 x 10 s 2.4 x 10 4 2.4 x 1(T 18 

plane EM wave (Cislunar distance) 



(b) 10 



-30 



3 x 10 10 Gaussian beam 



2.5 x 10 2 2.5 x 10 



-28 



(c) 10 



-30 



3 x 10 10 Gaussian beam 
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2.4 x 10 4 2.4 x 10" 



■26 



(d) 10~ 30 3 xl0 10 Gaussian beam 5 x 10 9 3.1 x 10 12 3.1 x 10~ 18 

(LISA dimension) 

V. CONCLUDING REMARKS 



1. For the relic GW's expected by the quintessential inflationary models, since a large amount of 
the energy of the GW's may be stored around the GHz band, using smaller EM systems (e. g., the 
microwave cavities or the Gaussian beam discussed in this paper) for detection purposes seems more 
plausible. In particular the Gaussian beams can be considered as new possible candidates. For the high- 
frequency relic GW with typical order of u g = 10 10 Hz, h = 10~ 30 in the models, under the condition 
of resonant response, the corresponding first-order perturbative photon flux passing through the region 
lCT 2 m would be expected to be 10 3 s _1 . This is the largest perturbative photon flux we have recently 
analyzed and estimated using the typical laboratory parameters. 

2. In our EM system, the perturbative effects produced by the + and x polarization states of the 
high-frequency relic GW have a different physical behavior. Especially, since the first-order tangential 
perturbative photon flux produced by the x polarization state of the relic GW is perpendicular to the 
background photon fluxes, it will be unique nonvanishing photon flux passing through the some special 
planes. Therefore any photon measured from such photon flux at the above special planes may be a 
signal of the EM perturbation produced by the GW, this property may be promising to improve further 
the EM response to the GW. 



3. As for the geometrical phase shift produced by the high-frequency relic GW, because of the 
excessive small amplitude orders of the relic GW, the phase shift is still below the requirement of the 
experimental observation. Thus the outlook for such schemes may be not promising unless there are 
stronger high-frequency relic GW's. But for the dimensions of the LISA (of course, in this case it needs 
a very strong microwave beam), it is possible to get an observable effect. 

The relic GW's are quite possible the few windows from which we can look back the early history 
of our universe, while the high-frequency relic GW's in the GHz band expected by the quintessential 
inflationary models, are possible to provide a new criterion to distinguish the quintessential inflationary 
and the ordinary inflationary models. As pointed out by Ostriker and Steinhardt [5], whatever the origin 
of quintessence, its dynamism could solve the thorny problem of fine-tuning our universe. If we could 
display the signal of the high-frequency relic GW's in the quintessential inflationary models through 
the EM response or other means, it would not only previde incontrovertible evidence of the GW and 
quintessence, but also give us the extraordinary opportunity to look back the early universe. Therefore, 
there is even a small chance that the signal of the high-frequency relic GW's is detectable, then it is 
worth pursuing. 

Finally, it should be pointed out that the relic GW's are the stochastic backgrounds, their signals 
have often larger uncertainty. However, it can be shown that for the relic GW's spectra expressed as 
Eq. (12), only one component (e.g. the form expressed as Eq. (15)) with special propagation direction, 
frequency, polarization and phase, can produce maximal perturbation to our EM system. Therefore, the 
high-frequency relic GW's would be instantaneously detectable at least. Moreover, if we can find a very 
subtle means in which the pure perturbative photon flux can be pumped out from our EM system, so that 
the perturbative and background photon flux can be completely separated, then the detection possibility 
of the relic GW's would be greatly increased. The more detail investigation will be discussed elsewhere. 
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APPENDIX A: THE CAVITY ELECTROMAGNETIC RESPONSE TO THE GRAVITATIONAL 

WAVES 



For the cavity electromagnetic response to the GW's, the best detecting state is the resonant response 
of the fundamental EM normal modes of the cavity to the GW's, since in this case it is possible to 
generate the maximal EM perturbation. Under the resonant states (whether the background EM fields 
stored inside the cavity is only a static field or both the static magnetic field and the normal modes), the 
displaying condition at the level of the quantum nondemolition measurement can be written as [10, 1 1] 

{hQ) 2 B 2 V > 

where Q is the quality factor of the cavity, V is its volume, B is the background static magnetic field. In 
order to satisfy the fundamental resonant condition, the estimated cavity dimensions should be compara- 
ble to the wavelength of the GW. However, the dimensions cannot be too big. This is because construct- 
ing a superconducting cavity with typical dimension I > 100m may be unrealistic under the present 
experimental conditions. The low-frequency nature of the usual astronomical GW's seem to greatly 
limit the perturbative effects in the cavity's fundamental EM normal modes. For the high-frequency 
GW's in GHz band, the corresponding resonant condition can be relaxed, but the cavity size cannot be 
excessive small, even if the condition I > \ g can be satisfied, since in this case the cavity cannot store 
enough EM energy to generate observable perturbation. If GW's detected by the cavity have excessive 
high-frequency, e. g., v g = v e > 10 9 Hz, we can see from Eq.(Al) that the requirements of other pa- 
rameters will be a big challenge. For instance, if one hopes to detect the high-frequency relic GW with 
h = 10~ 30 and v g = 10 9 Hz in the quintessential inflationary models, we need Q = 10 12 , B = 30T 
and V = 100m 3 at least, then the corresponding signal accumulation time will be r ps Q/co e w 10 3 s. 
If h = 10~ 30 , v g = 10 8 Hz, we need Q = 10 12 , B = 30T and V = 10m 3 (i. e., the typical dimension 
of the cavity will be I ~ 2.2m) at least, then r ps Q/co e ps 10 4 s. Increasing quality factor Q and 
using squeezed quantum states may be a promising direction [11]. For the former, the requirements of 
other parameters can be further relaxed; for the latter, the signal accumulation time could be decreased. 
Therefore, for the fundamental resonant response, the suitable size of the cavity may be the magnitude 
of meter, the corresponding resonant frequency band should be 10 8 Hz < v g < 10 9 Hz roughly (the 
region (4)-2 in Figure I). In order to detect the GW's of 10 8 Hz < v g < 10 9 Hz, the thermal noise must 
be KT < his g which corresponds to T < 10~ 3 fT (see Appendix D). 

1 1 



Moreover, the resonant schemes of difference-frequency suggested in Refs.[38, 39] can be used as 
EM detectors of the high-frequency GW's. The detector consists of two identical high-frequency cavities 
(e. g., two coupled spherical cavities discussed in resent paper [39]). When the GW frequency v g equals 
the frequency difference \v\ — v 2 \ of the two cavity modes (i. e., v g = \v\ — v 2 \, and v\, v 2 3> v g ), then 
detector can get maximal EM energy transfer. Following [39] we can learn that sensitivity of the EM 
detector would be expected to be 5h ~ 10~ 20 — 1CT 22 for the GW's of 10 3 — 10 4 Hz frequency band (the 
region (4)-l in Figure I). If this EM detector is advanced, it might to detect the GW's in GHz band. Ref. 
[40] reported an EM detecting scheme of the high-frequency GW's by the interaction between a GW 
and the polarization vector of an EM wave in repeated circuits of a closes loop. In the scheme because of 
linearly cumulative effect of the rotation of the polarization vector of the EM wave, expected sensitivity 
can reach up to Sh ~ 1(T 18 - 1(T 19 for the GW's of 10 8 - 10 9 Hz (the region (4)-2 in Figure I), and the 
above two scheme [39, 40] are all sensitive to the polarization of the incoming GW signal. Although 
the sensitivity of the above EM detectors is still below the requirements of the observable effect to the 
high-frequency relic GW's, the advanced EM detector schemes would be promising. 

APPENDIX B: THE DIMENSIONLESS AMPLITUDE h AND THE POWER SPECTRUM S h OF 
THE HIGH-FREQUENCY RELIC GRAVITATIONAL WAVES 

The relation between the logarithmic energy spectrum VLqw an d the power spectrum Sh of the relic 
GW given by Ref.[l] (see Eq.(A18) in Ref.[l]) is 

47T 2 

tt G w(v,r]o) = ^jp^S h (iy,r]o), (Bl) 

where H is the present value of the Hubble constant, i.e., H = 3.24 x 10 _18 s _1 . From Eq.(Bl), one 
finds [1] 

S h (v,m) ~ 8 x lO-^n GW (u, m )^^. (B2) 

In the peak region of the logarithmic energy spectrum of the relic GW in the quintessential inflation- 
ary models, Q GW « 5 x 1(T 6 [1]. Thus, for the relic GW of v = 10 9 Hz, we have S h « 4 x l(T 69 s; for 
the relic GW of v = 10 w Hz, one finds S h w 4 x l(T 72 s. 

For the continuous GW's, the dimensionless amplitude h can be estimated roughly as 

h^(S h Au)^ (B3) 

where Au g is corresponding bandwidth. According to estimation in Ref.[l] (see Figure.2 in Ref.[l] or 
see Figure. 1 in this paper), the bandwidth in the high-frequency peak region is about Av xs 10 11 — 



W 9 Hz « 10 11 Hz. Thus, from Eq.(B3), we have 

h « 1(T 29 ~ 1(T 30 (for the relic GW of 1/ = 10 9 Hz), 

h « 10 -30 ~ 10" 31 (for the relic GW of ^ = 10 10 Hz). (B4) 

The above results and orders of magnetude estimated in Ref.[l] are basically consistent. In this paper 
we have chosen h ~ 1CT 30 . Of course, these estimations are only approximate average effects. In fact, 
because of uncertainty of some relative cosmological parameters [1, 2] in certain regions, it is possible 
to cause small deviation to the above estimations. 

APPENDIX C: MINI-ASTRODYNAMICAL SPACE TEST OF RELATIVITY USING OPTICAL DE- 
VICES (MINI-ASTROD) 

Mini-ASTROD is a new cooperation project (China, Germany, etc.)[26]. The basic scheme of the 
Mini-ASTROD is to use two-way laser interferometric ranging and laser pulse ranging between the 
Mini-ASTROD spacecraft in solar-system and deep space laser stations on Earth to improve the pre- 
cision of solar-system dynamics, solar-system constants and ephemeris, to measure the relativistic 
gravity effects and to test the fundamental laws of space-time more precisely, to improve the mea- 
surement of the time rate of change of the gravitational constant, and to detect low-frequency GW's 
(~ 10~ 6 Hz - 10- s Hz). 

The follow-up scheme of the Mini-ASTROD is ASTROD[41], i.e., the Mini-ASTROD is a down- 
scaled version of the ASTROD. Both LISA and Mini-ASTROD are all space detection projects, but 
there are some differences for their study objectives, detecting frequency band of the GW's for the 
Mini-ASTROD will be moved to 10~ 6 Hz (see Fig.l). With optical method, Mini-ASTROD should 
achieve the same sensitivity as LISA [26]. Thus, Mini-ASTROD [26], ASTROD [41] and LISA have a 
certain complementarity. 

APPENDIX D: NOISE PROBLEMS 

The noise problems of the EM detecting systems have been extensively discussed and reviewd [1, 
10, 35, 38, 39], here we will give only a brief review to relevant problems of the EM detection of the 
high-frequency GW's (especially the EM response of the Gaussian beam). 

The thermal noise is one of the fundamental sourse of limitation of the detecting sensitivity [3 8, 39]. 
Unlike the usual mechanical detectors, the resonant frequencies of the EM systems to the high-frequency 



GW's in GHz band are often much higher than ones of usual environment noises (e.g., mechanical, 
seismic and others). Thus the EM detecting systems are easier to reduce or shield external EM noise 
(e.g., using Faraday cages) than closing out mechanical vibration from a detection system. For the 
EM response of the microwave cavities to the high-frequency GW's, the noise problem can be more 
conveniently treated considering the relevant quantum character [35]. For a superconducting cavity at 
a temperature T = T , if the background EM field is only a static magnetic or static electric field, the 
displaying condition can be given by Eq.(Al), while then the cavity vacuum contains thermal photons 
with an energy spectrum given by the Plank formula: 



where u u and v are energy density and the photon frequency respectively, while K is the Boltzmann 
constant. If the cavity is cooled down T = ImK, according to the Wien law, the energy density has a 
maximum at v m = 5.87 x 10 7 Hz (i.e., A m 3m, corresponding photon density is about lCT 8 cm~ 3 ). 
For the perturbative photons produced by the high-frequency GW of v g = 3 x 10 9 Hz under the resonant 
condition, we have v e = v g (i.e., A = 0.1m), which is higher than v m (i.e., v e « 30u m ). Therefor, the 
crucial parameter for the thermal noise is the selected frequency and not the total background photon 
number. Namely, in this case the thermal noise can be effectively suppressed as long as the detector can 
select the right frequency. 

For the EM response of the Gaussian beam in high-frequency region of v e = v g = 3 x 10 w Hz, 
because the frequencies of usual environment noise are much lower than u e , the requirements of sup- 
pressing such noise can be further relaxed. While for the possible external EM noise sources, using 
a Faraday cage would be very useful. Once the EM system (the Gaussian beam and the static mag- 
netic field) is isolated from the outside world by the Faraday cage, possible noise sources would be 
the remained thermal photons and self-background action. However, since the "random motion" of 
the remained photons and the specific distribution of the photon fluxes in the EM system (as we have 
discussed earlier), the influence of such noise to the highly "directional" propagated perturbative pho- 
ton fluxes would be effectively suppressed in the local regions. Therefore, the key parameters for the 
noise problems are the selected perturbative photon fluxes (e.g., and n',) passing through the special 
planes (e.g., planes <fi = n/2 and 3n/2) and not the all background photons. Moreover, low-temperature 
vacuum operation might effectively reduce the frequency of the remained thermal photons and avoid 
dielectric dissipation. If the frequency v m of the remained thermal photons is much lower than that of 
the perturbative photon fluxes, i.e., v m <C v e , then the above two kind of photons would be more easily 
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distinguished. 
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